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ABSTRACT 

We investigate the change in stellar magnetic topology across the fully-convective bound- 
ary and its effects on coronal properties. We consider both the magnitude of the open flux that 
influences angular momentum loss in the stellar wind and X-ray emission measure. We use 
reconstructed maps of the radial magnetic field at the stellar surface and the potential-field 
source surface method to extrapolate a 3D coronal magnetic field for a sample of early-to-mid 
M dwarfs. During the magnetic reconstruction process it is possible to force a solution towards 
field geometries that are symmetric or antisymmetric about the equator but we demonstrate 
that this has only a modest impact on the coronal tracers mentioned above. We find that the 
dipole component of the field, which governs the large-scale structure, becomes increasingly 
strong as the stellar mass decreases, while the magnitude of the open (wind-bearing) mag- 
netic flux is proportional to the magnitude of the reconstructed magnetic flux. By assuming a 
hydrostatic and isothermal corona we calculate X-ray emission measures (in magnitude and 
rotational modulation) for each star and, using observed stellar densities as a constraint, we 
reproduce the observed X-ray saturation atRo <0.l. We find that X-ray rotational modulation 
is not a good indicator of magnetic structure as it shows no trend with Rossby number but can 
be useful in discriminating between diff'erent assumptions on the field geometry. 
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1 INTRODUCTION 

Stellar coronal magnetic activity is predominantly investigated in 
the X-ray and radio bands. A linear relationship between the X-ray 
luminosity, Lx, and the radio luminosity, Lr, established by Guedel| 
|& Benz|(l993| ), pointed to a close connection between the hot coro- 
nal plasma responsible for the thermal X-ray emission and the free 
electrons causing the non-thermal radio emission, 

LogLx ^ LogLR + 15.5. (1) 

While this linear relationship holds for several classes of ac- 
tive, main sequence star, between types F and early-M, and is also 
verified for T Tauri stars and solar flares, the assumption that the 
Lx - Lr relationship would continue into the realm of the very low- 
mass stars proved to be short lived; the detection of radio emission 
from an M9 dwarf, LP944-20, defied this relation by almost four 
orders of magnitude (^Berger et al.||200r) . Since the discovery of 
radio emission from LP944-20, many more observations have been 
carried out on low-mass stars to reveal that they also show the same 
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radio phenomenon fBerger"2002; "Berger et al."2005", 'Bu rgasser &| 
|Putman 2005 , Hallinan et al. 2006 ; McLean et al. 201 1 ). 

Radio observations are important diagnostic tools as they can 
help determine the magnetic field configuration ('Halli nan et al.| 
[2008 ) and the nature of plasma emitting regions. X-ray observa- 
tions on the other hand help provide a deeper understanding of 
the generation of the magnetic field. Originating principally in the 
magnetically confined plas ma of the hot stellar corona - at tem- 
peratures greater than 10^^ ('G udel|20'04| - X-ray emissions prove 
useful for determining stellar parameters, such as the extent of the 
stellar corona, which also aids in establishing the structure and evo- 
lution of the magnetic field. X-ray emission has shown to be preva- 
lent throughout the main sequence, with the maximum value of 
X-ray luminosity observed for each spectral type decreasing with 
bolometric luminosity, (e.g. |Barbera et al.|1993j ) - for spectral type 
down to early-M. 

K and M stars are the most numerous in the Galaxy, and cou- 
pled with their high levels of X-ray emission, they are prime candi- 
dates for investigating the dependence of X-ray emission on physi- 
cal parameters such as mass, rotation rate, and more recently, mag- 
netic topology. So far, the X-ray luminosity has been shown to cor- 
relate well with either rotational velocity or Rossby number (the ra- 



2 P. Lang et al 



tio of the stellar rotation period, P, to the convective turnover time, 
Tc) (e.g. Mangeney & Pra derie||1984l [MarilH & Catalano|[T984l 
[Schmitt et al.||1985| [Fleming et al.||19 89 ); in genera l, Lx/Lhoi in - 
creases and then saturates (Lx/Ltoi ~ 10~^; ( Delfosse et al.|l998] )) 
with increasing rotation rate. This behaviour is attributed to coronal 
saturation ( Vi lhu & Walter|1987[[Stauffer et al.|1994| ); however, it 
is not clear whether saturation is a reflection of the dynamo itself 
or the total filling of the stellar surface with active regions ( |Vilhu| 
[T984J . 

For a small sample of stars later than M5, the Lx/Lboi ratio was also 
shown to decline once again at the lowest Rossby numbers within 
the sa mple stars (e.g.[Gol ub 1983, Bookbinder 1985; Rosne r et"ar] 
p^85l|Jemes et al.|2011| ) and [Berger| ( 2006) found that the value 
drops to Lx/Ltoi ~ 10~^^ at spectral type M7. This indicates that 
outwith flaring, the emission from the thermal plasma cannot ex- 
ceed a certain fraction of the total stellar flux. This super-saturation 
has been attributed to negative dynamo feedback ( Kitchati nov et al.| 
1 1994) 1, lack of coverage of active regions ( Stgpien et al.|200l| ), or 
centrifugal stripping of the corona ( [Jardine 2004[ ). 

It is important to note that there is a change in internal struc- 
ture between the early-M dwarfs (dM) and mid-M dwarfs. The 
former have a solar-like internal structure: a turbulent, electrically 
conducting convective envelope, surrounding a radiative core. The 
interface between these two zones i.e. the tachocline, is the site 
where amplification of the magnetic field is believed to take place 
([Spiegel & Zahn' 1992V The relative size of the radiative core dra- 
matically drops with decreasing temperature for early M dwarfs, 
and at spectral type M4 O.SSM©), the stars are fully convec- 
tive (e.g. [Chabrier & Baraff'e|1997] ). This change in internal struc- 
ture coincides with a change in magnetic topology ( D onati et al.| 
[20081 IMorin et al.||2008l |20T0l hereafter D08, M08 and Ml 0, re- 
spectively). In the higher mass stars, M ^ O.45M0, the field con- 
figuration is more complex, i.e. a non-axisymmetric poloidal field 
with a strong toroidal component, in comparison to the nature of 
the field structure at spectral types later than M4, which is mainly 
axisymmetric and poloidal. 

As for magnetic activity, it is crucial to look at coronal prop- 
erties such as emission measure, especially for fast rotators, to ex- 
plore whether saturation and super- saturation occur at a fixed rota- 
tion period or Rossby number. This could potentially help in deter- 
mining the physical mechanism that is responsible for the satura- 
tion, i.e. an eff'ect of the dynamo efficiency, or a consequence of the 
fast rotation rates exhibited by these stars. 

Stellar models (e.g. |Gilliland|1986[ [Kim & Demarqu e|1996| ) 
suggest that the convective turnover time is longer in lower mass 
stars; if one considers the Sun and an M-dwarf with an equivalent 
rotation period, then the M dwarf has a smaller Rossby number and 
yet is more active in terms of Lx/Ltoi- Although coronal saturation 
occurs at Lx/Lhoi ~ for spectral type G,K and M, the rotation 
period at which it sets in is larger for the lower masses. Therefore, 
this results in a Rossby number of 0.1 being the point where coro- 
nal saturation sets in ( [Patten & Simon' 1996';Pizzolat o et al.|2003[ 
[Jeff^ries e t al. 201 li However, Jeffries et al. (2011) show that the 
super- saturation is more clear as a function of rotation period within 
each spectral type, rather than Rossby number, with the eff'ect oc- 
curring at periods of P < 0.3 days for K-dwarfs and P < 0.2days 
for M dwarfs. 

Magnetic activity and rotation are well correlated ( [Stauff'erj 
[etal.[1997| and both are strongly dependent on stellar age for main- 
sequence solar-like stars ( [Skumanich|1972| ), 



where Qq is stellar angular velocity. During their lifetime, rota- 
tional evolution of stars can be governed by disc braking, pre-main- 
sequence contraction i.e. spin-up, and/or magnetic winds. From ob- 
servations of young open clusters (e.g jlrwin et al.|2009j , it is clear 
that rapidly rotating dM stars are a common occurrence; however, 
in older open clusters, this number is reduced. [Scholz et aL] ( [2011| ) 
find a mass-dependent exponential rotational braking law 
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where towards lower masses the spin-down timescale r increases, 
with T ^ 0.5Gyr for 0.3Mo and r > IGyr for O.IM©. This indicates 
that at an age of 600Myr stars of 0.3Mo are almost exclusively fast 
rotators. 

As well as X-ray emission, another coronal tracer of magnetic 
activity is radio emission. The relation that correlates the X-ray 
luminosity with the radio luminosty (eqn. [T]), suggests that ther- 
mal X-ray emission, assumed to be due to hot coronal plasma, 
and non-thermal radio emission, possibly generated by the elec- 
tron cyclotron maser instability ( |Melrose & Dulk|1982] ) or gyrosyn- 
chrotron, are both consequences of the magnetic field. Throughout 
the subtypes of dM stars, the radio luminosity remains approxi- 
mately constant, while Lx tracks the bolometric luminosity (i.e. 
Lx/Ltoi ~ 10~^). The sharp deviation from the Lx - Lr relation, 
first noted by jBerger et al. l('200r), does not present itself until be- 
yond spectral subtype M7 (Berger et al. 2010 ), where the relation 
evolves from Lr/Lx ~ 10~ 



to ^ 10 This would indicate 
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that the deviation must not directly relate to the transition to full 
convection. Along with this increase in radio luminosity, it is still 
unclear as to why the emission is present on these stars during one 
set of observations and then absent during the next (e.g. Berger| 
[et al.|2010| ). 

In order to investigate the role of field topology on the coronal 
structure and emission properties of M dwarfs, we use the observed 
surface magnetic field maps of their coronae. From this we can 
predict the X-ray emission measure. 



2 SAMPLE SELECTION 

Spectropolarimetry along with multi-line extraction methods (e.g. 
least- squares deconvolution), has the power to detect Zeeman sig- 
natures in the polarisation spectra of low-mass stars. With a sig- 
nificant improvement in instrumentation and imaging capabilities, 
we are now able to recover information on the distribution of large- 
scale magnetic fields on the surface of cool active stars (^Donati] 
[et al. [12003] ). This has allowed for a much more in-depth investi- 
gation of dynamo processes in low-mass stars. Using rotationally 
modulated circular polarisation signatures, the photospheric mag- 
netic field can be reconstructed from a series of polarised spectra, 
through Zeeman-Doppler Imaging (ZDI). The magnetic field is de- 
composed into its toroidal and poloidal components and described 
as a set of spherical harmonic coefficients ( [Donati et al.|2006[ ). 

The first low-mass star on which this was attempted was the 
M4 dwarf V374 Peg, which exhibited a very strong, large-scale, 
mainly axisymmetric poloidal field ([Donati et al.|2006[ >, which is 
stable on timescales of ^ lyr ( [Morin et al.|2008| >. This result lead 
D08, M08 and MIO to investigate dM stars further: a survey was 
conducted using data from ESPaDOnS and NARVAL on 23 active, 
main- sequence, dM stars spanning a broad range of masses and ro- 
tation rates. The mass range embodied dM stars located both above 
and below the fully convective boundary and the rotation periods 
were derived using ZDI, which provides a constraint on the rotation 
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period as well as the projected equatorial velocity (v sin /). Rossby 
number are from |Donati eTaLl ( [2008] ); rMorin et al.| ( [2008] ) and were 
computed from empirical suited to the stellar mass from |Kiragal 
|& Stepien |([2007]). 

In this paper, we concentrate on the M-dwarf sample with 
masses ranging from 0.25Mo to 0.75Mo (DOS, M08), where 
0.75Mo refers to a very young star GJ182, and compare results to 
stellar parameters such as mass and Rossby number as these can 
provide the best insight into dynamo action and efficiency. 



3 FIELD GEOMETRY 

3.1 Potential-Field Extrapolation 

To calculate the magnetic field above the stellar surface, a potential- 
field source surface (PFSS) approach is applied ( [Altschuler &| 
|Newkirk|| 19691 . By assuming that the coronal magnetic field is 
current-free (V x B = 0), the magnetic field can be defined as the 
gradient of the scalar potential, B = - Vj/^. Since the magnetic field 
is divergence free (V • 5 = 0), this results in Laplace's Equation: 
V^if/ = 0. The solution to Laplace's equation in spherical harmon- 
ics allows for the determination of the coefficients of the associ- 
ated Legendre polynomials P/^ and, therefore, the coronal mag- 
netic field components. 
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where aim and bim are the amplitudes of the spherical harmonics, / 
is the spherical harmonic degree and m is the order or 'azimuthal 
number' . 

To extrapolate the 3D coronal field and determine the ampli- 
tude of the spherical harmonics, aim and bim, we apply two bound- 
ary conditions. The upper condition is that at some distance from 
the star i.e. the Source Surface, /?^^( |Schatten et al.|1969| ), the field 
opens and is purely radial {Bq - B^p - 0). The amplitude of the 
spherical harmonics aim and bim can then be written as 



Cllm 
bim 



1 
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where Rss is chosen to be 2.5P* {R^ is the stellar radius), based 
on the value for the solar source surface and for simulations of 
V374Peg ( |Vidotto et al.|20Tl] ). 

The lower boundary condition completes the extrapolation 
by imposing the observed radial field to determine the values of 
bim for all values of / (up to 31 - for model resolution purposes) 
and m considered. The code used to extrapolate the field is a 
modified version of the global diff'usion model developed by |van| 
|Ballegooijen'etar| ( [7998l ). 



3.2 Stellar Surface Maps 

Fig.( |l|2| and[3]l show the reconstructed maps of the surface radial 
magnetic field, viewed from our line of sight with respect to the 



stellar inclination, /(°) (see Table [TJ, for the early-to-mid M dwarf 
sample. The ZDI maps show significant structure at the stellar sur- 
face (a and b). This degree of structure decreases with height until 
on the largest scales the extrapolated fields are well described by 
tilted dipoles (c). Since DOS, M08 and MIO report strong compo- 
nents of the poloidal field at low Rossby numbers, we investigate 
the polar field strength of the dipole component. 

Fig. |4] shows the polar field strength of the dipole component 
to increase with decreasing Rossby number. Using this combined 
with Table [T] we find that the most noticeable transition occurs in 
the stars with mass < 0.40Mo, where the field strength increases 
from < 200G to ^ 400G. This behaviour follows the results 
obtained by DOS, MOS & MIO and indicates that the large-scale 
field has a stronger dipole component for fully convective stars. 



3.3 Surface Flux and Open Flux 

In order to investigate the coronal structure of our sample of stars, 
we determine the variations with Rossby number of both the total 
magnetic flux through the stellar surface and also the amount of 
that flux that is open. The open flux is of particular interest because 
it channels the stellar wind. Open field lines are therefore not only 
dark in X-rays (and hence contribute to any rotational modulation 
of the X-ray luminosity) but they are also responsible for the angu- 
lar momentum loss by the stellar wind. 

The ratio of open flux to surface flux is given by 



Rl j \Br(R,,e,(p)\dQ ' 



(8) 



S urface 



where dCl is the solid angle. 

To analyse the geometry of the field, which is described by 
the overall combination of / and m modes, we compare the pre- 
dicted total open flux of each star in the sample to that obtained by 
only single value modes. By choosing only single modes, we can 
simplify equation JSl to 



_ (2/+l)(^)^^ 
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(9) 



and compare the predicted open-to-observed surface flux, calcu- 
lated through the combination of / and m modes for each of our 
stellar sample, to 3 single modes: a dipole (/=i), a quadrupole {1-2) 
and an octupole {1-3). 

From Fig. [5] (top) we can deduce that although (near) fully- 
convective objects have increasingly large-scale, low multipole 
field topologies (DOS, MOS & MIO), demonstrated by the extrap- 
olated 3D coronal fields (Fig[T][2]and[3]) and the strong polar field 
strengths of the dipole component (Fig. |4]), we cannot ignore the 
higher multipole components when modelling magnetic activity, or 
even the magnetic field itself. For example, if the dipole compo- 
nent, which has the most open flux associated with it, was assumed 
to be the only significant magnetic field component, the amount of 
open flux could be overestimated by at least an order of magnitude. 
This would have c onsequences for, say a Weber-Davis model fWe-| 
ber&Davis|l967| ) of angular momentum loss, /, due to the stellar 
wind. 
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which could overestimate the angular momentum loss by 2 orders 
of magnitude. We note here that the magnitude of the open flux 
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(a) GJ182 



(a) DT Vir 



(a) DS Leo 



(a) GJ 49 



(a) OT Ser 




Figure 1. (a) and (b): The reconstructed surface radial maps shown at longitudes 180° apart to give the best viewing angle of both magnetic poles, (c): The 
extrapolated field, using the PFSS method with Rss = 2.5R*. For each star, colours are scaled to the maximum and minimum values of the radial magnetic 
field component: blue represents negative flux and red positive flux. 
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(a) CE Boo 




(a) AD Leo 



(a) EQ Peg A 






(a) EV Lac 



(b) 






(a) YZ CMi 



(b) 




Figure 2. (a) and (b): The reconstructed surface radial maps shown at longitudes 180° apart to give the best viewing angle of both magnetic poles, (c): The 
extrapolated field, using the PFSS method with Rss = 2.5R*. For each star, colours are scaled to the maximum and minimum values of the radial magnetic 
field component: blue represents negative flux and red positive flux, (cont.) 
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(a) V374 Peg 






(a) EQ Peg B 



(b) 




Figure 3. (a) and (b): The reconstructed surface radial maps shown at longitudes 180° apart to give the best viewing angle of both magnetic poles, (c): The 
extrapolated field, using the PFSS method with Rss = 2.5R*. For each star, colours are scaled to the maximum and minimum values of the radial magnetic 
field component: blue represents negative flux and red positive flux, (cont.) 
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tempting to extrapolate from this result that, according to a Weber- 
Davis model ( [Weber & Davis |p"967| , the increased flux on these 
stars should cause a high rate of angular momentum loss, due to 
the stellar wind (Eq.[To]), we lack an appropriate stellar wind model 
that would enable us to calculate this loss and make any definite 
conclusions. Furthermore, without stellar ages for our sample, it 
is difficult to determine whether these stars follow an age-rotation 
law e.g. ([2| or a mass dependent rotational braking law (Eqn[3]), i.e. 
have high rotation rates because they are young. 

[Vidotto et al.| ( |20 lfl modelled the angular momentum loss of 
the M4 dwarf V3 74 Peg through its stellar wind, and in conjunction 
with its rapid rotation rate, V374 Peg exhibits a high angular mo- 
mentum loss, suggesting it is still young. However, for comparison 
purposes, an approximate stellar age for V374 Peg has not yet been 
determined. 



Figure 4. Polar field strength of the aligned dipole component i.e. (m=0, 
1=1) as a function of Rossby number where asterisks represent M > OAMq 
diamond symbols represent M < OAMq. Although no distinct transition 
is observed across the coronal saturation boundary (Ro = 0.1), it is clear 
that stars approaching full-convection (diamonds) have strong polar dipole 
components. Note that the outlier in this plot is M0.5 dwarf DT Vir (GJ 
494A) at Ro =0.0092. 

is heavily dependent on the chosen value for the source surface as 

Rss =^ ^Open 0. 

Fig. [5] (bottom) demonstrates that, as with the polar field 
strength (Fig. |4|, the amount of surface and therefore open flux 
increases with inverse Rossby number, albeit with a more notable 
flattening at the fully convective boundary. This indicates that the 
lower mass stars (diamonds) exhibit a significant fraction of open 
flux [lO^'^Mx] and supports the idea that the large-scale dipole com- 
ponent is having an eflfect on the magnetic activity. Although it is 



4 MODELLING THE EMISSION 
4.1 The X-ray Emission Model 

Once the global magnetic field is determined, the gas pressure and 
density structure of the stellar corona can be estimated. We assume 
the plasma is hydrostatic and isothermal (T = 2x 10^ K; Giampapa] 
|et al.| ( p^96| ) and that the gas pressure at the stellar surface is pro- 
portional to the magnetic pressure (po oc B^). The pressure there- 
fore varies along each field line according to 

r s:Bds 

P^Poe^ 1^1 , (11) 

as described by^ Jardine et al.| ( [2002| and [Gregory et al.| ( |2006| ). 
Expanding the component of gravity along the field line (g • B), we 
have 
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Figure 5. Top: The magnitude of the predicted open flux (where asterisks 
represent M > 0.4Mo and diamond symbols represent M < 0.4Mo) as a 
function of the observed surface flux is influenced by the higher multipole 
field components. Although the open flux tracks the surface flux, which has 
a large dipole component {1=1), the overall combination of / and m modes 
is closer to that of an octupole Bottom: The open flux (shown in 

blue) and the surface flux (shown in black) as a function of Rossby number 
shows both to increase with lower Rossby number. This follows the same 
trend as the polar field strength of the dipole component shown in Fig.|4]but 
there is a more pronounced flattening at the onset of coronal saturation i.e. 
Ro <Q.\. Note that again the outlier in this plot is M0.5 dwarf DT Vir (GJ 
494 A). The stars shown coloured are discussed in detail in section|5] 
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where /c is a constant of proportionality and Br.Bg^B^ are the ra- 
dial, meridional and azimuthal components of the magnetic field, 
respectively. Finally, the ratios of centrifugal-thermal energy (0^) 
and gravitational- thermal energy (0g) are given by 



[ ksT j 



(13) 



(14) 



where Q is the stellar rotation rate, kg is the Boltzmann constant, G 
is the gravitational constant and is the electron mass. 

To ensure that only regions of the closed stellar corona con- 
tribute towards the emission measure, the gas pressure along open 



field lines is taken to be zero. In addition to this, if there is any 
over-pressure along the designated closed loops (/? ^ 1 .0), then the 
pressure at that grid point is also set to zero. Here p is the ratio of 
gas pressure {p - lUekT) to magnetic pressure {pB - B^/2fi). 

The X-ray emission measure is then determined assuming the 
gas is optically thin, i.e. 



EM{r) : 



(15) 



and we consider both its magnitude (Fig.|6| and rotational modula- 
tion i.e. variation in emission measure (Fig.|7]). 

The high levels of rotational modulation indicate that even on 
stars where the magnetic field is mainly dipolar, i.e. low Rossby 
number, there is still enough structure present in the magnetic field 
to give significant rotational modulation. Of course, we must also 
consider the large variation in stellar inclination angle and angle 
between the magnetic axis and rotation axis, i.e. /3m (see Table [TJ, 
when analysing the modulation to determine whether it is a good 
indicator of magnetic structure. From Fig. ( |l|2| and[3]l and Table 1, 
we can see that the magnetic axis and rotation axis are mis-aligned 
by /5m > 10° for the vast majority of the sample and this, coupled 
with the stellar inclination angles, could account for the large vari- 
ation between rotational modulation values within the sample. 

The free parameters in this model are the temperature (fixed 
at r = 2 X 10^^), the source surface (fixed at Rss = 2.5/?*) and 
the constant of proportionality k relating the base gas to magnetic 
pressure. Changing k essentially scales the global plasma pressure 
up or down and we use observations of coronal density to fix k. 
[Ness et al] ( |2002| |2004| ) estimates coronal densities on the dM 
stars to be within the range of 10^ - lO^^cm"^. Fig. |6] (upper) 
shows the maximum and minimum density occurring in the closed 
corona of each sample star. There is a great spread in the range 
of maximum-to-minimum densities calculated for each star, with 
some extending over as much as 4 orders of magnitude within the 
closed corona. However, we note that these density ranges represent 
diflferent (r, 6, 0) values for each star. The average coronal density 
calculated in our model for each star, which can be found in Table[T] 
are estimated to lie within the range 10^ - lO^^cm"^, representing 
early-dM to mid-dM respectively. We note the former would appear 
to encompass densities comparable to solar values. 



5 EXTRAPOLATION ASSUMPTIONS 

In the reconstruction of the surface magnetic field it is possible to 
push the solution towards field geometries that are either symmet- 
ric (e.g. a quadrupole) or antisymmetric (e.g. a dipole) about the 
centre of the star (see for example Fig[8]). The magnetic modelling 
is described in detail in DOS. The principle of ZDI requires a way 
to choose between the many field configurations that fit the data 
equally well. The imaging process is based on the principles of 
maximum entropy image reconstruction and in the absence of any 
prior information on the field topology, an entropy form that selects 
the solution with the lowest energy content and the lowest degree 
modes, is generally used. Should we have a good physical reason to 
think that the field is predominantly symmetric or antisymmetric, 
we could impose such a constraint with a modified entropy. In com- 
parison to the unconstrained solution, magnetic energy is higher 
(and entropy lower) for both the symmetric and anti- symmetric 
cases because the magnetic field is on smaller spatial scales. We 
remind the reader that the antisymmetric solution is not only dipo- 
lar. The ZDI solution with the lowest energy (weighted by the or- 
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Figure 6. Top: The maximum and minimum coronal density occuring on 
each the stars in the sample (where the long- dashed line represents M < 
0.4Mo and the short-dashed line represents M > 0.4Mo. Note the spread 
in the density values calculated for each star can extend over many orders 
of magnitude but represent different values for (r, 6, 0). Bottom: Emission 
Measure as a function of Rossby number (where asterisks represent M > 
OAMq and diamond symbols represent M < 0.4Mo) . This demonstrates 
the trend observed as the boundary of coronal saturation is approached i.e. 
Ro < 0.1. Note the short dashes are the partly convective stars and the long 
dashes are the fully convective stars. 



der of the spherical harmonics) that fits the data is selected. When 
one forces antisymmetry, the number of degree of freedoms in the 
model is decreased, and hence, if the unconstrained solution was 
not purely antisymmetric, the antisymmetric solution has a lower 
entropy/higher energy. We investigate the eff'ect of this choice on 
coronal extrapolations (and quantities derived hereafter). In partic- 
ular, we determine if the X-ray emission measure is sensitive to the 
nature of the constraint since comparison with the observed emis- 
sion measure may then help to discriminate between choices in the 
way that for T Tauri stars, the location of accretion spots does ( Do- 
Inati et al.|20TT| ). 

We explore the assumptions made for surface field images of 4 
stars: MO, DS Leo; M1.5, OT Ser; M3 AD Leo; and M4 V374 Peg, 
which encompass high, mid and low inclination angles, as well as 
a range of stellar masses. Figure |8] shows the coronal extrapolation 
for the 3 cases: anti- symmetric, symmetric and unconstrained, for 
V374 Peg. We note here a few remarks on the reconstruction of 
the individual stars. For AD Leo, which has a low inclination an- 
gle, 70°, a predominantly dipolar field is reconstructed with ZDL 



Figure 7. The rotational modulation of the X-ray emission as a function of 
Rossby number (where asterisks represent M > 0.4Mo diamond symbols 
represent M < OAMq) depicts the variation in emission measure during a 
full rotation of the star. Even at low Rossby number, there is still enough 
structure present in the magnetic field to give significant rotational modula- 
tion. 



However, for the field in the invisible hemisphere instead of a clear 
pole of strong magnetic field, it is covered with field of average 
intensity with polarity opposite to that of the visible pole. As a re- 
sult the field lines from the visible pole connect on the whole hid- 
den hemisphere whereas for a pole-like model they would connect 
on the well-defined opposite pole. By imposing an antisymmetric 
field this eff'ect is corrected for: the resulting ZDI map features a 
strong polar spot in the invisible hemisphere with opposite polarity 
with respect to the visible pole. For DS Leo, the inclination angle 
is intermediate and the map does not feature any strong spot on 
the partly hidden hemisphere. It is likely in this case, as for AD 
Leo, that we lack information and do not know what is on the less- 
visible hemisphere. Imposing antisymmetry indeed results in spots 
of high magnetic field on this hemisphere, and also modifies more 
importantly the map in the visible hemisphere. For a star with a 
well constrained dataset and high inclination angle like V374 Peg, 
the map is already well constrained and features an overall dipo- 
lar structure with a strong polar spot on the hemisphere orientated 
away from the observer. 

Whether the unconstrained or antisymmetric solution is forced 
we find a similar result in terms of the predicted open flux (Fig. [9]) 
and the magnitude of X-ray emission measure. Imposing a symmet- 
ric field, one artificially imposes smaller spatial scales e.g. /^m = 2, 
i.e. the lowest degree mode available becomes the quadrupole, 
which shifts magnetic energy from order / to (at least) l+l and au- 
tomatically decreases the amount of open flux. This can be consid- 
ered as a test on a star for which a good dataset is available, e.g. 
V374 Peg. We note that even in the cases of forcing a solution (ei- 
ther symmetric or antisymmetric), the relation between the open 
flux and surface flux is still followed. 

The eff'ect of imposing diff'erent solutions has very little or 
no eff'ect on the X-ray Emission Measures or coronal densities 
calculated. We do however note a change in the X-ray rotational 
modulation between solutions (Fig.[To]). An observation of the X- 
ray rotational modulation (such as that described in |Hussain et aL| 
( |2004| |2005 ), may discriminate between fields that are predomi- 
nantly symmetric or antisymmetric. As we discussed in Section |4] 
this depends on the inclination of the magnetic pole and the inch- 
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Figure 9. The variation in the open flux as a function of the surface flux 
when difl'erent configurations of the magnetic field are assumed. The solu- 
tions are represented as follows: antisymmetric, crosses; symmetric, square. 
The unconstrained solutions are represented by asterisks if M > 0.4Moand 
diamond symbols if M < OAMq). Artificially imposing a symmetric or anti- 
symmetric field increases the reconstructed magnetic energy and hence both 
the observed surface flux and predicted open flux. The symmetric solution 
(squares) has the lower open flux, which is to be expected since imposing a 
symmetric field shifts magnetic energy from order / to (at least) l+l. 



(b) Symmetric 




(c) Unconstrained 

Figure 8. Coronal extrapolation of the (a) antisymmetric, (b) symmetric, 
and (c) unconstrained ZDI solutions of the coronal extrapolation for V374 
Peg. Note that only the closed field lines are shown. 



nation of the star itself. For some stars e.g. V374 Peg, this is not so 
useful, whereas for other like AD Leo there is a large variation in 
the rotational modulation making it a good candidate. 
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Figure 10. The eff'ect of imposing diff'erent magnetic field symmetries 
on the X-ray rotational modulation. The solutions are represented as fol- 
lows; antisymmetric, crosses; symmetric, square. The unconstrained solu- 
tions are represented by asterisks if M > 0.4Moand diamond symbols if 
M < OAMq). Although we have shown that rotational modulation is not a 
good indicator of field structure when considering a large sample of stars 
with diff'erent inclination, etc, it could be useful when considering a single 
object, e.g. checking the symmetry and providing independent confirmation 
of the predominant mode. 



6 SUMMARY AND CONCLUSIONS 

We have used reconstructed maps of the radial magnetic field at the 
stellar surface for a sample of early-to-mid M dwarfs to extrapolate 
their 3D coronal magnetic field (using the PFSS method). We have 
investigated the topology of the large-scale magnetic field at the 
stellar surface and the structure of the extrapolated 3D corona. By 
assuming a hydrostatic, isothermal corona, we have modelled the 
density structure within the corona and hence determined the X-ray 
emission measure. 
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We have focussed in particular on variations with Rossby 
number. We find the following: 

1. As the Rossby number decreases, the polar field strength 
of the dipole component of the field increases and then appears 
to saturate. Stars with low Rossby numbers have strong, mainly 
dipolar fields. 

2. A similar variation with Rossby number is seen in both the 
total (unsigned) surface magnetic flux and the flux of open field 
(which can carry the stellar wind). The role of the topology of the 
large-scale field is apparent when we calculate the magnitude of 
the open flux. This is significantly less than would be predicted if 
all of the surface magnetic flux were contained in a purely dipole 
field. The contribution of the higher multipoles therefore reduces 
the open flux and may also significantly influence the angular mo- 
mentum loss rate, which for a Weber-Davies model scales as the 
square of the open flux. Both the strength and also the topology of 
the large-scale field are therefore important in angular momentum 
loss. 

3. As is observed, a rise and then saturation of the X-ray emis- 
sion measure with decreasing Rossby number is also found. The 
stellar coronae are compact, with most of the emission originat- 
ing from regions below approximately 1.5 stellar radii. Our sample 
contains a large range of both the inclinations of the stellar rotation 
axis and also the tilt of the magnetic axis. As a result, there is a 
large spread in the values of the rotational modulation of the X-ray 
emission and no clear trend with Rossby number can be detected. 

For low-mass stars, the observed variation in X-ray emission 
with Rossby number results naturally from the observed variation 
in the surface magnetic field. We note that we choose the parame- 
ter K that scales the surface pressure in such a way that we repro- 
duce typical X-ray fluxes and we do not model the ionisation frac- 
tion sometimes invoked in the atmospheres of later spectral types. 
While there is a range of magnetic topologies within our sample, 
the spread of values for the rotational modulation of the X-ray 
emission is too great for it to be a useful indicator of the field struc- 
ture. Although not a good indicator of field structure when consid- 
ering a large sample of stars with diff'erent inclination, it could be 
useful when considering a single object. For example, checking the 
symmetry and providing independent confirmation of the predomi- 
nant mode (e.g. dipole versus quadrupole), or observing a magnetic 
cycle with a change from predominantly quadrupolar to predomi- 
nantly dipolar as on the Sun { [Sanderson et al.|2(j03] ). 

We find that both the strength of the field and its geometry, 
however, aff'ect the magnetic flux that is open (wind-bearing) and 
which therefore allows the star to lose mass and angular momen- 
tum. The magnitude of this open flux is significantly reduced by 
departures from a purely dipolar field. This suggests that simple 
scalings for angular momentum losses based on dipolar field ge- 
ometries may not be sufficient to explain the angular momentum 
evolution of low mass stars. The high values of open flux for stars 
with the lowest Rossby numbers may indicate that they have sig- 
nificant angular momentum loss rates. 
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Table 1. Data for stellar sample of early-to-mid M dwarfs. Mass, radius, rotation period, inclination and, where available. By, the average large-scale magnetic 
field derived from spectropolarmetric measurements are provided by Donati et al. ( 2008 1; Morin et al. ( 2008 1. Bi, is the average magnteic field (i.e. small -i- 
large-scale field) derived from unpolarised spectroscopy, supplied by a) Reiners & Basri ( 2009 1, b) Saar (1996), c) Reiners & Basri (2007), d) Johns-KruU &] 
|Vale nti ( 2000). Rossby number are fromlDonati et al.ll2008); Morin et al. ( 2008) and were computed from empirical Tc suited to the stellar mass from Kirag^ 
|& S tepien ( 2007). fi, the estimated angle between the rotation and magnetic axis, are from this paper, along with the predicted values for emission measure 
(both magnitude and rotational modulation) and coronal density. 



Star 


Sp Type 


Mass (Mq) 


Radius (Rq) 


P (days) 


Ro(10-2) 


in 


/3m C) 


By (kG) 


Bi (kG) 


LogEM (cm ^) 


Rot Mod 


Log^e (cm ^) 


GJ 182 


M0.5 


0.75 


0.82 


4.35 


7.44 


60 


41.1 


0.172 


2.5^ 


50.33 


12.18 


8.62 


DT Vir 


M0.5 


0.59 


0.53 


2.85 


9.2 


60 


83.6 


0.149 


3.0^ 


50.91 


2.82 


9.26 


DS Leo 


MO 


0.58 


0.52 


14.0 


43.8 


60 


41.1 


0.087 




48.36 


9.93 


7.98 


GJ 49 


M1.5 
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0.51 


18.6 


56.4 


45 


10.9 


0.027 




46.83 


18.02 


7.05 


OTSer 


M1.5 
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0.49 


3.40 


9.7 


45 
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0.123 




50.66 


48.83 


9.22 


CEBoo 


M2.5 
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14.7 
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7.4 


0.10 
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49.47 


2.17 


8.52 


AD Leo 


M3 


0.42 


0.38 


2.3399 


4.7 


20 


4.5 


0.19 


2.9^ 


50.22 


3.27 


8.95 


EQ Peg A 


M3.5 


0.39 


0.35 


1.061 


2.0 


60 


25.9 






51.42 


51.02 


9.66 


EV Lac 


M3.5 


0.32 


0.30 


4.3715 


6.8 


60 


45.8 
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52.02 


12.05 
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YZ CMi 
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> 3.9' 


50.21 


9.89 


10.27 


V374 Peg 
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0.44565 
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70 


9.3 
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26.49 


10.26 


EQ Peg B 


M4.5 


0.25 


0.25 


0.404 


0.5 


60 


6.5 






51.11 


14.05 


9.64 



jen A. A., Collier Cameron A., Jardine M., Donati J.-R, 2004, in 
A. K. Dupree & A. O. Benz ed., Stars as Suns : Activity, Evo- 
lution and Planets Vol. 219 of lAU Symposium, Using Chan- 
dra/LETG Spectra to Probe Stellar Coronae. p. 301 

Hussain G. A. J., Brickhouse N. S., Dupree A. K., Jardine M., 
van Ballegooijen A., Collier Cameron A., Donati J.-R, Pavata 
R, 2005, in R Ravata, G. A. J. Hussain, & B. Battrick ed., 13th 
Cambridge Workshop on Cool Stars, Stellar Systems and the Sun 
Vol. 560 of ESA Special Publication, Inferring coronal structure 
using X-ray spectra: a Chandra study of AB Dor. p. 665 

Irwin J., Bouvier J., Hebb L., Scholz A., Stassun K., 2009, in 
NOAO Proposal ID #2009B-0144 Rotation periods of solar- 
mass stars in NGC 2547 and NGC 2516. pp U4-+ 

Jardine M., 2004, A&A, 414, L5 

Jardine M., Wood K., Collier Cameron A., Donati J.-R, Mackay 
D. H., 2002, MNRAS, 336, 1364 

Jeffries R. D., Jackson R. J., Briggs K. R., Evans R A., Pye J. R, 
2011, MNRAS, 411, 2099 

Johns-Krull C. M., Valenti J. A., 2000, in R. Pallavicini, 
G. Micela, & S. Sciortino ed.. Stellar Clusters and Associations: 
Convection, Rotation, and Dynamos Vol. 198 of Astronomical 
Society of the Pacific Conference Series, Measurements of stel- 
lar magnetic fields, pp 371-+ 

Kim Y.-C, Demarque R, 1996, ApJ, 457, 340 

Kiraga M., Stepien K., 2007, ACTAA, 57, 149 

Kitchatinov L. L., Ruediger G., Kueker M., 1994, A&A, 292, 125 

Mangeney A., Praderie R, 1984, A&A, 130, 143 

Marilli E., Catalano S., 1984, A&A, 133, 57 

McLean M., Berger E., Reiners A., 2011, ArXiv e-prints 

Melrose D. B., Dulk G. A., 1982, ApJ, 259, 844 

Morin J., Donati J.-R, Porveille T, Delfosse X., Dobler W, Petit 
P., Jardine M. M., Collier Cameron A., Albert L., Manset N., 
Dintrans B., Chabrier G., Valenti J. A., 2008, MNRAS, 384, 77 

Morin J., Donati J.-R, Petit R, Delfosse X., Porveille T., Albert L., 
Auriere M., Cabanac R., Dintrans B., Pares R., Gastine T., Jar- 
dine M. M., Lignieres P., Paletou R, Ramirez Velez J. C, Theado 
S., 2008, MNRAS, 390, 567 

Morin J., Donati J.-R, Petit P., Delfosse X., Porveille T., Jardine 
M. M., 2010, MNRAS, 407, 2269 



Ness J.-U., Gudel M., Schmitt J. H. M. M., Audard M., Telleschi 
A., 2004, A&A, 427, 667 

Ness J.-U., Schmitt J. H. M. M., Burwitz V, Mewe R., Raassen 
A. J. J., van der Meer R. L. J., Predehl R, Brinkman A. C, 2002, 
A&A, 394,911 

Patten B. M., Simon T., 1996, ApJS, 106, 489 

Pizzolato N., Maggio A., Micela G., Sciortino S., Ventura P., 
2003, in A. Brown, G. M. Harper, & T. R. Ayres ed.. The Puture 
of Cool-Star Astrophysics: 12th Cambridge Workshop on Cool 
Stars, Stellar Systems, and the Sun Vol. 12, Activity-rotation Re- 
lationship: Interpretation of an X-ray Derived Rossby Number, 
pp 887-892 

Reiners A., Basri G., 2007, ApJ, 656, 1121 

Reiners A., Basri G., 2009, A&A, 496, 787 

Rosner R., Golub L., Vaiana G. S., 1985, ARA&A, 23, 413 

Saar S. H., 1996, in K. G. Strassmeier & J. L. Linsky ed.. Stellar 
Surface Structure Vol. 176 of lAU Symposium, Recent magnetic 
fields measurements of stellar, pp 237-+ 

Sanderson T. R., Appourchaux T., Hoeksema J. T., Harvey K. L., 
2003, Journal of Geophysical Research (Space Physics), 108, 
1035 

Schatten K. H., Wilcox J. M., Ness N. R, 1969, Sol. Phys., 6, 442 
Schmitt J. H. M. M., Golub L., Harnden Jr. R R., Maxson C. W, 

Rosner R., Vaiana G. S., 1985, ApJ, 290, 307 
Scholz A., Irwin J., Bouvier J., Sipocz B. M., Hodgkin S., 

Eislofifel J., 2011, MNRAS, 413, 2595 
Skumanich A., 1972, ApJ, 171, 565 
Spiegel E. A., Zahn J.-R, 1992, A&A, 265, 106 
Staufifer J. R., Hartmann L. W, Prosser C. R, Randich S., Bal- 

achandran S., Patten B. M., Simon T., Giampapa M., 1997, ApJ, 

479, 776 

Staufifer J. R., Liebert J., Giampapa M., Macintosh B., Reid N., 

Hamilton D., 1994, AJ, 108, 160 
St?pien K., Schmitt J. H. M. M., Voges W, 2001, A&A, 370, 157 
van Ballegooijen A. A., Cartledge N. R, Priest E. R., 1998, ApJ, 

501, 866 

Vidotto A. A., Jardine M., Opher M., Donati J. P., Gombosi T. L, 

2011, MNRAS, 412, 351 
Vilhu O., 1984, in E. Rolfe ed., Pourth European lUE Conference 



12 P. Lang et al 



Vol. 218 of ESA Special Publication, Magnetic activity in lower 

main sequence stars, pp 239-242 
Vilhu O., Walter R M., 1987, ApJ, 321, 958 
Weber E. J., Davis Jr. L., 1967, ApJ, 148, 217 

This paper has been typeset from a T^X/ OTgX file prepared by the 
author. 



